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TECHNICAL NOTE 3^36 


THE ACCURACY OF THE SUBSTITUTE-STRINGER APPROACH FOR 
DETERMINING THE BENDING FREQUENCIES OF 
MULTISIRINGER BOX BEAMS 
By William W. Davenport 


SUMMARY 


The accuracy of the substitute-stringer approach for including the 
effects of shear lag in the calculation of the transverse modes and fre- 
quencies of multistringer box beams is investigated- Box beams , the covers 
of which consist of normal-stress -carrying stringers on sheets which carry 
not only shear but also normal stress , are analyzed exactly. Frequencies 
of beams with various numbers of stringers, obtained by means of this 
exact analysis, serve to determine the possible accuracy of the frequencies 
obtained by the substitute-stringer approach. 


INTRODUCTION 


The use of a substitute-stringer approach for including the effects 
of shear lag in the vibrational analysis of built-up box beams was Inves- 
tigated in reference 1. Various thin-walled rectangular tubes were Ideal- 
ized to substitute-stringer structures and the frequencies of the first 
few bending modes of the idealized structures were compared with those 
of the original beams. The results indicated where the substitute stringer 
should be located in order to take into account accurately the effect of 
shear lag on the frequency; for a fairly wide variety of tube proportions, 
locating the substitute stringer midway between the web and the centroid 
of the half -cover yielded frequencies accurate within 1 or 2 percent for 
higher as well as lower modes of vibration. However, because the rec- 
tangular tubes employed in the investigation had constant wall thickness 
there was no definite - Inf ormation given regarding the accuracy of fre- 
quencies that may be obtained for a beam which has discrete flanges and 
stringers in the covers and different web and cover-sheet thickness. 

The purpose of this report is to examine the accuracy of the 
substitute-stringer approach when applied to a beam more nearly repre- 
sentative of actual built-up box beams. Frequencies of various realistic 
multistringer beams of the type shown in figure l(a) are determined by 




2 


KACA TN 5636 


means of the substitute-stringer approach as given in reference 1. These 
frequencies are then compared •with the frequencies obtained by means of 
an exact vibrational analysis of a multistringer beam. Conclusions are 
thereby made with regard to the proper location of the substitute stringer. 


SYMBOLS 


A 

Ap 

a l 

*0 

A p 

A S 


a 

b 

b 0 = a 
b P 

b r+l 


b S 

b C 

Bi('l) 


C 

Ci (cl) 


Di 


(0 


total cross-sectional area of stringers in a half -cover of 
multi stringer box beam 

cross-sectional area of flange of substitute-stringer structure 
cross-sectional area of a substitute stringer 
cross-sectional area of flange of multistringer box beam 
cross-sectional area of pth stringer of multistringer box beam 
effective shear-carrying area 

half -depth of beam 
half -width of beam 

width of pth panel 

half -width of middle panel when IT is even, nonzero constant 
when ...K is odd 

distance between web and adjacent substitute stringer 
distance between web and centroid of area of half -cover 

parameter defined by equation (AI 7 ) or equation (kko) 
constant 

parameter defined by equation (Al 8 ) or equation (A4l) 
parameter defined by equation (A25) or equation (A^ 2 ) 
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E 


modulus of elasticity 




G 


SVihq. 


Gt C L 2 ( L + 0 N S (r+l),q) 


shear modulus of elasticity (taken equal to E/2.65 herein) 


H = Sl 

t c a 


bending moment of inertia 


ki/ 4 ’ - + 6 oj) + <J*) ! 


S B 


frequency coefficient , a> 





coefficient of shear rigidity. - \E- 

lVgas 


k i = " k B 2k S £ 


L 

N 

Ni 


(P) 


half-length of beam 

number of stringers in a cover 

parameter defined by equations (A34) and (A35) or by equa- 
tions (Alf3) and (Abh ) 


P q. “ 


m r 


GL 


Ri (cl) = Fq 2 (^) 2 + BS 


iq. 


r = |(N + 0 N ) 
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- 

Sl (p) 

determinant given "by equation (A3Q) 




cover-sheet thickness 

-■ 

. - 

t c' 

effective cover-sheet thickness for normal stress 

y 

tw 

web thickness 



T 

maximum kinetic energy 



U 

maximum strain energy 




longitudinal displacement of a point on 

the pth panel 

— 

w 

vertical displacement of cross section 

of beam 

— 

X 

longitudinal coordinate 


- 

y P 

chordwise coordinate for pth panel 

- — 

.. V l 


Fourier series coefficients 


w 

integers 



B = b( l 




Hc l 2Gt c L 




8 13 

Kronecker delta (0 when i 3} 1 when 

i - J) 


e n = li - 

<-n N ] 



M (p) „ 

Lagrangian multiplier 



n 

mass of "beam per unit length 



*l (p) 

parameter defined by equation (A8) 


* 

01 

natural frequency 


u 

oie 

natural frequency of multi stringer beam 

in figure 3 
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THE MULTISTRINGER BOX BEAM AND THE SUBSTITUTE-STRINGER STRUCTURES 


A doubly symmetrical multistringer box beam of the type considered 
in the present paper is shown in figure l(a). In order that the covers 
of this structure may behave realistically with regard to shear-lag 
effects , the cover sheets are permitted to carry not only shear but also 
normal stress. The following simplifying assumptions are made: 

(1) The flanges and stringers carry only normal stress. 

(2) The webs carry only Shear (the bending resistance of the web 
is included in the flanges). 

(3) Longitudinal inertia (rotary inertia) is neglected. 

(4) The cross sections maintain their shape. 

Assumptions (2) and (3) are generally known to be good for reasonably 
shallow beams such as those considered in this report; assumption (4) 
is good if a normal amount of bulkhead stiffness is present. 

In the appendix an exact vibration analysis, similar to those of 
references 1 and 2, is carried out for a multistringer beam of the type 
shown in figure 1(a). The Rayleigh-Ritz energy procedure is used in 
conjunction with appropriate Fourier series and Lagrangian multipliers 
to obtain frequency equations for a box beam with any number of stringers. 
Although the analysis in the appendix allows unequal spacing of the 
stringers, only the case of equally spaced stringers is considered 
hereinafter . 

The substitute-stringer structure for the beam of figure l(a) is 
shown in figure l(b). The flange-web combination of this structure is 
the same as that of its prototype; the covers, however, consist of sub- 
stitute stringers which carry only normal stress and sheets which carry 
only shear. The magnitude of the shear-lag effect in this structure 
depends on the location of the substitute stringers. A suggested stringer 
location given in reference 1 is 


b s = 0.5bc (l) 

where bg is the distance between a flange and the adjacent substitute 
stringer and bg is the distance between a web and the centroid of 
normal-stress -carrying area of the half -cover. 

A vibration analysis of the substitute -stringer structure is carried 
out in appendix B of reference 1. 
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COMPARISON OF FREQUENCIES AS OBTAINED BY THE MULTISTRINGER 
SOLUTION AND THE SUBSTITUTE-STRINGER APPROACH 


In order to Investigate the possible accuracy of the substitute- 
stringer approach for determining the frequencies of multistringer box 
beams , comparisons are made between the frequencies of four multistringer 
box beams as obtained by the analysis given in the appendix and the corre 
sponding frequencies of their substitute-stringer structures as obtained 
by the analysis given in reference 1. The first three symmetrical trans- 
verse modes of each beam with free-free end conditions are considered. 


The four beams have cross sections as shown in figure 2 and are 
identical except for the manner in which the total stringer area is 
distributed. Proportions common to all four beams are 



trj 

-T = 1.25 
t C 


~ = 1-89 
at w 


~~ = 0.75 

btc 


where L, b, a, t^, and t^ are, respectively, the half-length, the 

half -width, the half -depth, the web thickness, and the cover-sheet thick- 
ness; the total cross-sectional area of the stringers in a half-cover is 
given by A and the cross-sectional area of a flange is given by Aq. 

In addition to the preceding properties, the four beams are alike in that 
stringers of each beam are equally spaced across the covers. The beams 
thus differ only in the number of stringerB on a cover; the cases con- 
sidered are for N = 1, 4, 7 } and 00 . 


In figure 2 the various relations needed to obtain the parameters 
of the multistringer solution for the four beams are given. When N Is 
Infinite the stringers are smeared out to a "stringer sheet" for which 
the effective normal-stress-carrying thickness tQ 1 is different from 

the actual cover-sheet thickness tg. 

The cross section of the substitute-stringer idealization Is also 
shown In figure 2. In accordance with the Idealizing procedure of the 
substitute -stringer approach, the following relations are used to deter- 
mine the parameters needed in the substitute-stringer solution of refer- 
ence 1: 


Ap = Aq Al = A + btQ 

where Ap and Aj, are, respectively, the cross-sectional area of a 
flange and stringer of the substitute-stringer beam. It should be noted 
that the only geometric property of the substitute-stringer structure 
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which may he different for each of the four teams is the value of bg. 

For the substitute-stringer location defined by equation (l), bg varies 
by virtue of the variation of bg. The values of bg for the four beams 
are: 


For 1 stringer. 


For If- stringers. 


bg = 0.714b 


bg = 0.543b 


For 7 stringers. 


b c = 0.531b 


For infinity stringers, 


b c = 0.500b 

The frequency coefficients k B = to 

substitute-stringer solution where bg = 0.5bg and the multistringer 

solution for each of the four beams are given in table I. In addition 
to these results, the frequency coefficients obtained by a solution 
which includes transverse shear but not shear lag have been included. 
This latter set has been presented in order to demonstrate the influence 
of shear lag and was obtained frcm equation (15) of reference 3 with the 
rotary inertia parameter equal to zero. 

The influence of stringer location on the frequency co of the 
substitute-stringer structure for each of the four beams i6 presented 
graphically in figure 3* The frequency is given in the form of its 

relative error ^ - 1 when compared with the exact frequency co e of 

the multistringer structure’, the curves indicate how the errors vary 
with substitute-stringer location. 



obtained from both the 


DISCUSSION OF RESULTS 


The results given in table I show good agreement between the fre- 
quency coefficients as obtained by the substitute-stringer approach with 
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bg/bg =0.5 add those obtained by the multistringer solution. However, 
as may be noted, the frequency coefficients for the substitute structure 
with bg/bg =0.5 are slightly high for all but the third modes of the 

cases for N = 1 and N = °°. Reference to figure 3 indicates that, for 
all but the case of N = <», the maximum errors for the modes considered 
would be reduced by using a value of bg/bg = 0.55> for example, for the 

case of N = 7, the maximum error for this new value would be less than 
1 percent . 

The results of table I indicate that the stringer area distribution 
has little influence on the frequency; however, the assumption should 
not be made that the influence of shear lag is negligible. As may be 
seen in table I, the reduction in frequency due to the inclusion of shear- 
lag effects is 4.8, 16.7, and 19.1 percent, respectively, for the first 
three symmetrical modes of the case where N = 7* 


CONCnmiRG REMARKS 


The numerical results of the present paper indicate that the 
substitute-stringer approach can yield accurate frequencies for multi- 
stringer box beams. Review of the numerical results of MCA Technical 
Note 3158 together with those of the present paper suggests that a 
value of O.55 for bg/bg (ratio of the distance between web and adjacent 

substitute stringer to distance between web and centroid of the area of 
the half-cover) defines a slightly more appropriate substitute-stringer 
location for built-up box beams than bg/bg = 0.50, the value suggested 
in MCA Technical Note 3158. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 6, 1958. 
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APPENDIX 

VIBRATION SOLUTION OF A MULTISTRINGIR BOX BEAM 


t 


* 


The simplifying assumptions for the multistringer beam are given 
in the body of the paper and a cross section consisting of two cases is 
shown in figure If. Case I applies when N is even and case II when 
N is odd, where N is the number of stringers in a cover. In ad d i tion 
to the dimensions and coordinates defined by figure if, the length of the 
beam is assigned the value 2L, and x is defined as the longit udina l 
coordinate with its origin at the midpoint of the beam. It should be 
noted that b r+ p, the distance between the rth stringer and the center 
line for case I, is half the width of the middle panel. 

For a transverse mode of vibration let w(x) be the amplitude of 
a vertical displacement of a cross section and let Up(x,yp) be the 

amplitude of a longitudinal displacement of a point on the pth panel. 
Then the maximum strain energy is 


- g c - -<•*«. ( i *H 


dx + 


20t , 


£ (i • CC ** 


dx + 


(s - 


(Al) 


where 


a 2 


N 


(A2) 


and 


5 (r+l),p - 0 
6 (r+l),(r+l) = 1 


(p i r + 1) (A3) 
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The maximum kinetic, energy when the influence of longitudinal inertia 
is neglected is , 



(A4) 


where p is the mass per unit length of the structure which the multi- 
stringer heam represents. It should he noted that equations (Al) to (Ah ) 
hold for N = 0. 

In the following sections appropriate infinite trigonometric series 
are assumed to represent the displacements and then the Rayleigh-Ritz 
energy procedure is applied. 


Symmetrical Transverse Modes of a Free -Free Multistringer Beam 

Appropriate trigonometric series for the symmetrical transverse 
modes of a free-free heam are 


V*,y P > - £ £ e n» (p)eln 1? «» ^ 

m=l,3,5 n=0,l,2 


2L hp 

(p * 1,2, . . . (r + 1)) 


w(x) = C + 


cos 


m-1,3,5 


rojtx 

2L 


(A5) 

(A6) 


The choice of these particular series was guided hy the orthogonality 
required for simplification of the expressions in the strain-energy 
equation; C in equation (A 6) is introduced in order to allow w(±L) 
to he unrestricted. 


Because of the geometry of the structure the following restraining 
relations must hold: 

up(x/bp) - u p+1 (x,0) =0 (p = 1,2, (AT) 

Substituting equation (A 5 ) into equation (A7)> multiplying hy sin 
where i = 1,3>5> • ■ • , and integrating from -L to L yields 


T. L (p) (-U n - a ln (p+1 3 - «*i (p) - 0 

-r\ i n — — 1 


(p “ . • . r)(i = • • •) 


(AS) 
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In accordance with the Rayleigh-Ritz procedure the expression 

U - T - XI 51 (A9) 

P=l,2 m=l,3,5 

is held stationary with respect to the coefficients of the assumed series. 
The A's in expression (A9) are Lagrangian multipliers introduced to 
maintain the condition stated in (A8). The expressions for U and T 
of equation (A9) are given by equations (Al) and (Ak), respectively, and 
the displacements are in turn given by equations (A5) and (A 6). Thus, 

differentiating expression (A9) with respect to the aj^^'s, ci's, 
and C independently and setting the respective results equal to zero 
yields the following equations: 


Kij (L Wj (l) + 2Ri (L ' ) 51 + Y HltY 1 ' Pl^i (1) (- 1 ) J = 0 

n=0, 1,2 L 


(l) x 2a Tr /ijt'i 


(1 ),_,J = 


(i — 1,3,5, • • • )(j — 0, A,2, . . .) 


(A10) 


(! - 0 #(r + l),g)%3 (q> «lJ ( ' l) + 2Kl (,l) X a ln (< ^ > - 

n=0 $ 1 j 2 

Pq|Ai (<l) (-l) J (l - S(r+l),q) - = 0 

(i = 1,3,5, • • • ) ( j = 0,1,2, . . . ) (q = 2,3, • • • r + l)(r = l) (All) 

kiCi + X a iM (l) - 2^^) (-1)^0 = G 

' ' n-0,1,2 

(i = 1,3,5, • • •) (A12) 

00 rci-I 

k B 2 C+ X Cmf^C-l)"^ - =0 (A13) 

_ m=l,3,5 

where 


Kij (tl) = P q 2 (|^) 2 (l + 5 0J ) + (inf 


(Mk&) 
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2 Eb n c 

p * -1& 


R i 


(4) = F„ 2 m z + 


^ V2; 


B5 


Id 


^g-l^q 


Gt G L 2 (i + e N s (r+I) ^) 


H = ^i 
t c a 


2 _ 






El 


*8=-^ 

GAgL^ 

Ag = 4aty 


3q 


- b d 


2Gt C L 


= 


{%)" - %V 


Dividing equation (AIO) by and summing over j 

(1). 


gives , on solving for y a^_» 

J=0A;2 


3 = 0 , 1,2 


1J (1) - B t (: %H (1) - $ 


(i - 1,3,5, • • 


(Al4b) 

(AIM) 

(Al4d) 

(Alke) 

(AIM) 

(Al4g) 

(AIM) 

(AIM) 

(AIM) 

0,1,2, « • • 
0 (A15) 
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Similar ly, for equation (All): 


£ H 3 M = ^i (4) (l - S( r+1)q )- c 1 ( ' 1 V q >>i (4 ' 1) 

j=0 ,1,2 ^ 


■where 


(i = 1,3,5; 


(-1 ) l 


)(q = 2,3, • • . r + 1) (A16) 


Bi 


(q) = 3=0, 1,2 


1 + 2Ri 


(4) 


3=0,!,2 K ^ (q) 


(a = 1,2, 


r) (AIT) 


and 


Ci 


(4) _ 


3=0, 1 , 2 K i,1 


nr 


(! - 0 #( r + i ), q ) + 2 % 


(4) 




3*071,2 K ± j 

(q = 1,2, . . . r + l) (Al8) 


In equation (AIT) the definition of BjL^ r+ ^ has "been neglected because 

when q = r + 1 in equation (Al6) the entire term containing 
drops out. 

By substituting equations (AI 5 ) and (Al6) into equations (AlO) 
and (All) the following equations may be obtained: 


n (1) — ^ 

Kij U; 


L(|)H(f) [fRi (l) Ci (l) - lj+ P 1 A i (l) j(- 1 ) 5 - 2R i (l) B i (l ^| 


(i = 1,3,5, - • -)(3 = 0,1,2, . . .) (A19) 
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K i.J 


tqT 1 


Pq4 


i (ci- 1 )[^ i (q.)c 1 (q.).^ 



(i = 1 , 3 , 5, • • 00 = 0,1,2, . . . )(q. = 2,3, . . . r) (A20) 

Equations (A19) and (A20), when multiplied by (-1)3 and summed over 
j = 0,1,2, . . yield 




(i = 1 , 3 , 5 . • • • ) (A2l) 


and 


S~ (-l^a*^ = + $q[ci (cl) + 7*^ 

j=0,l,2 “ J 

(i = 1,3,5, • • •)(<!“ 2,3, . . . r) (A22) 

In order to obtain equations in the form of (A2l) and (A22) it should be 
noted that 




-Bi 


(q) 


(q = 1,2,3, • • • f) (A23) 


and 


£ 1 £ M)L 



+ Di 


O) 


(q = 1,2,3, . • . r) (A24) 
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where 


2Ri 


M 


= 


5~ \ 

3=071,2 


(-D 


3' 


J=o,i,g K i3 


TiT 


/ \ go 

i + 2Rl w x 


3=0, 1,2 KlJ (<i) 

(q. = 1;2,5, . . - r) (A25) 

Substituting equation (A15) into equation (A12) results in 

[*1 - + - ^ 2 % 2 (£)c-i>~ £ 


(i = 1,3,5, • • •) (A26) 

and substituting equations (A2l) and (Al6) into equation (a 8), when 
p = 1, gives 


_B. ^ ^ ( ^ aN \ n fi^ \ C i 

1 \L ) \2 ) p L 


Ci (L) + Di (l) + ^ c^ 2) ; 


*i 


(D - 


^B^ 2 \ (2) (l- Hr+lh2 ) =0 


b l 

(i * 1,3,5, • • •) 

Substituting equations (Al6) and (A22) into equation (A8) yields 


(A2T) 


^(pyp-n + 


(P) . -n (P) . ^P+l „ (P+l) 


+ D. 


Cl 


(p) 


b^T B i (P+L)> i (P+l) (l " »(r+l),(p+l)) " 0 


(i = 1,3,5, . . - )(p = 2,3, . . . r) (A2B) 


By Cramer’s rule, equations (A26) to (A28) may be solved for C i/Pq 
in terms of C/|3q; if in the resulting expression the determinants in 



the numerator and denominator are each expanded hy the first column and if in the resulting 
denominator the cofactor of -B^ ^ ^ a -t- fl0 expanded by the first column, the following 

equation results: 


OV-p \( -1 \ ^ _£L a. (l) 


where 




•e 

- 2H(#)V ir 


* ) 

CfcW + Djbl . OjtP* 1 ) 


• ■ • 1 * 

• 



... 

■ 

■ ■ 1 


ftjO**) +i^0^«) + ^« I !»t») ... . 

• ■ i • * • 

f 

• 

. 

• ... « 

... Cl !r^) + Il (r-«) + i = l 0l (r- 1 ) 
■ive 

• 

• 

• 

... .. ^fr-U^^r-D^ 


• 

• 

■% W 

4 w, 4 w t to,H 



! ! 

[f - l.e.J, . . . r) WO) 

1 

1 

1 

i 

1 



k 

I .Mi/ 


I. . i 
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By expanding the determinant of (A30) hy the first col umn the following 
recurrence relation may be obtained: 

8i (p) - 0i.( p) + d/ p) + ^±i Cl (p+1 ^ - ^3: j^ ± Cp+iJ| S Si Cp+2 ) 

(p = 1,2,3, ... r - 1) (A31) 

where 

Si (r> ' [ci (r) + Di (r) ♦ ^±1 Cj.^' 

L_ °r _ 

and 

Bi (r+1) = 1 

Now, by dividing the numerator and denominator of the right-hand member 
of equation (A29) by and then using equation (A3l), the following 

equation results: 


i-l 



(1 = 1,3,5, • • •) (A32) 

Repetition of this procedure ultimately yields the continued fraction 
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permits equation (A53) to be written as 

„ = 2 C 


ki - 


(¥) 2 * 


(i = 1,3,5, • • •) (A36) 


■tc 


Hi 


(1) 


Substituting equation (A36) into equation (A13) yields the frequency 
equation 


k 2 

2 2 , a- 

1 4. \ 


1 


% < 

L + Kg \ — 

/ fa 

Ht\ 2 

/mrt\ 2 Hi 

. W ° 



L- 

m=l,3,5 

2 J 

^ ( 1 ) 

v _ 



- 0 (A3T) 


(P) 


The following closed forms for the summations involved in B m , 

Cn/ P ^, and Dm.^^ may be obtained by methods similar to those given in 
reference k: 

(p = 1,2, . . . r + l) (A38) 


0=4 t ,2 2^P p tanh 


L=±l 


A 


(pi 


5=°A,2 2(f-)Pp sinh (f-)Pp 


(p = 1,2, . . . r + 1) (A39) 


By use of equations (A38) and (A39), equations (AI7), (Al8), and (A25) 
may now be written as 


Bm 


(P) _ 


sech 


/mtrtV 
V 2/ 


(f)Pp (f)Pp + % 


(pH 


(p = 1,2, . . . r) (Ak0) 


^m 


(p) 


(p) 


(L - 9 #(r + l),p)(f) P P taDh (f) P P + ** 

(p = 1,2, . . . r + l) (Akl) 
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(p) _ V 


tanh W\? x 


2 (^) P P (f)Pp tanh 00P p + bJ 


(p = 1,2,3, • • • r) (Ak2) 

From equations (A3 4), (A35), (A40), (A4l), and (A42) the following expres- 
sions for and may be obtained: 


lim (p) = ili ^P p tanh (SL)Pp + I k (p) + 

^ ^2/ * \ 2 / tanh (SljP, 


(p = 1,2, . . . r)(b 0 » a ) (A43) 


(r+ 1 ) 


0 n) O ?) 1 ^ 1 t£mh (f) Pr+i + ^ (r+1) 


(r ^ 1) (A44) 

The rate of convergence of the series of frequency equations given 
in equation (A 37 ) is increased by substracting the expression 


4 2 
2% kg 


=1^,5 (M) kjn 


and adding the equivalent closed-form expression 


gftan^kgkg y 

B V k B*S " V 


The resulting equation is 


2 3 3 

k B tan kg kg + 2kg kg 


” =1 ' 5 ' 5 v* | % (1) 


- ^(fj 


= 0 (A45) 
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where is determined, "by equations (Alt-3) and (Aid)-) . It should be 

noted that equation (Alj-5) will hold, for N stringers where N / 0. 


For the case when N = 0, handling equations (A12), (A13)^ and (A10) 
when = 0 in a manner similar to that used for the case of 

N stringers yields equation (At 5 ) except that 



- £ ?i(f ) *i(f) 



(m) 


The values Pp and are given in equations (Al4b ) and (AlAc ), 

respectively; the value of b-j_ is, of course, b. 


In order that equation (Alt-5 ) together with (At-6) may apply to the 
case for N = °°, P^ may be redefined by 



Eb^t^’ 

GL 2 t c 


(M-7) 


Antisymmetrical Transverse Modes of a Free-Free Multistringer Beam 

Appropriate trigonometric series for the antisymmetrical transverse 
modes of a free -free beam are 


-QQ_- , 93 , 

u p (x,y p ) = 2_ Z_ 

m=0,2,4 n=0, 1,2 


w(x) = Cx + 

m=2 ,b,6 


^ cos cob 

rutyp 


2L 



(p = 1,2, . . 

. r + l) 

(Alt8) 

„ rnicx 

c m sin 


(Alt-9) 


As in the case for the symmetrical modes the choice of the particular 
series was guided by the orthogonality required for simplification of 
the expressions in the strain-energy equation; the term Gx was added 
in order to allow sufficient freedom of w(±L). 
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Treatment similar to that accorded the, case of the symmetrical modes 
yields the equation 


k n 


■B 


1 - k B kg cot k B kg + 2k B kg 




m= 2 , 4,6 ^2 ^ ^(1) . 


= 0 (A50) 


where kg, kg, and are given in equations (Alkf), (Alkg), and (Alkj), 

respectively, and is given by equations (Ak 3 ) and (Akk). As 

before , this equation holds for N stringers , where N / 0 . When N = 0 , 

equation (A50) may be used, where is defined by (Ak6); for N = », 

_L 1 _ __ _ — JL^ -> ^ . JXV VN-fi-4 /N-P P *1 /V>t V> _ 

respectively; the value of b B is, of course, b. 

In order that equation (Ak 5 ) together with (Ak6) may apply to the 
case for N = °°, P-^ may be redefined by 


Eb^tp 1 

p = ^ 

1 GL 2 t n 


(AkT) 


Antisymmetrical Transverse Modes of a Free-Free Multistringer Beam 

Appropriate trigonometric series for the antiBymmetrical transverse 
modes of a free -free beam are 


-QQ_- , CP , 

u p (x,y p ) = 2_ , Z_ 

m=0,2,k n=0, 1,2 


w(x) = Cx + 

m=2,k,6 


^ COS COS 

nityp 


2 L 



(p = 1,2, . . 

. r + l) 

(Ak8) 

_ nmx 

c m sin 


(Ak 9 ) 


As in the case for the symmetrical modes the choice of the particular 
series was guided by the orthogonality required for simplification of 
the expressions in the strain-energy equation; the term Gx was added 
in order to allow sufficient freedom of w(±L) . 



NACA TN 3656 


23 


REFERENCES 


1. Davenport, William W. , and Kruszewski, Edwin T. : A Substitute-Stringer 

Approach for Including Shear-Lag Effects in Box -Beam Vibrations. 

NACA TN 3158, 1954. 

2. Budiansky , Bernard, and Kruszewski, Edwin T. : Transverse Vibrations 

of Hollow Thin-Walled Cylindrical Beams. NACA Rep. 1129, 1953* 
(Supersedes NACA TN 2682.) 

3. Kruszewski, Edwin T.: Effect of Transverse Shear and Rotary Inertia 

on the Natural Frequency of a Uniform Beam. NACA TN 1909, 1949* 

4. Adams, Edwin P., and Hippisley, R. L. : Smithsonian Mathematical 

Formulae and Tables of Elliptic Functions. Second reprint, 
Smithsonian Misc. Coll., vol. 74} no. 1, 1947} p- 129* 



2k 


MCA TN 5636 


TABLE I 

EFFECT OF SPRINGER DISTRIBUTION ON FRE$JENCJy 


N = 4 


N = 00 


Solution 

k B 

1st symm. 
mode 

2d symm. 
mode 

3d symm. 
mode 

Elementary 

+ 

transverse shear 

5-42 

24-5 

48.3 

Substitute- 
stringer, b s /b c =0.5 

5.16 

20.1 

_i 

37-8 

Multi stringer 

5.07 

19.7 

38.1 

Elementary 

+ 

transverse shear 

5-42 

24.5 

48.3 

Substitute- 
stringer, bg/b(j = 0-5 

5.22 

20.8 

39-1 

Multi stringer 

5-15 

20.2 

38.8 

Elementary 

+ 

transverse shear 

5-42 

24.5 

48.3 

Substitute- 
stringer, tig/bc = 0-5 

5-22 

20.9 

39-2 

Multistringer 

3.16 

20.4 

39.1 

Elementary 

+ 

transverse shear 

5-42 

24.5 

48.3 

Substitute- 
stringer, "bg/t>c = 0.5 

5-23 

21.0 

39-5 

Multistringer 

5.19 

20.8 

39-9 
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Figure 3.- Effect of substitute-stringer location on the accuracy of the 
substitute -stringer approach for box beams with different stringer 
distributions. The labels IS, 2S, and 3S denote symmetrical modes. 
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Figure 4 .- Cross section of doubly symmetrical multistringer box beam 
for either even or odd number of stringers in a cover. 
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